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Abstract 

We present a method, based on loop equations, to compute recursively, all the terms 
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energy of a statistical physics model on a discretized torus. 



E-mail: eynard@spht.saclay.cea.fr 



1 Introduction 



Random matrix models [TP|, ^P], [FT], [37|] have a wide range of applications in 

mathematics and physics where they constitute a major field of activity. They are 
involved in condensed matter physics (quantum chaos [29, 34jl, localization, crystal 
growths P0|1,...etc), statistical physics |1^,[1^,[14],[28|] (on a 2d fluctuating surface, also 
called 2d euclidean quantum gravity, linked to conformal field theory), high energy 
physics (string theory ||, quantum gravity |19|, pj, |27| , QCD f44j],...), and they are 
very important in mathematics too: (they seem to be linked to the Riemann conjec- 
ture p6| , |38[| ), they are important in combinatorics, and provide a wide class of integrable 
systems 0,0, 

In the 80's, random matrix models were introduced as a toy model of O-dimensional 
string theory and quantum gravity [17, H], 14]. By fine-tuning the parameters of the 
model, one can reach many multi-critical points, which are in relationship with the 



(p, q) minimal conformal models |2U, [L6| . The 1-hermitian matrix models is limited 
to (p, 2) conformal minimal models, whereas the 2-matrix model can represent any 
(p, q) conformal minimal model |16| . The 2-matrix model is thus more general than 
the 1-matrix model, and it has been the source of considerable interest in the past few 
months [[||2g,|g,[lj],§,|,|l|^,[46 . 



The free energy of the 1-matrix model was conjectured [f£3j 19, |, |3| and also rigor- 
ously proven |22j to have a 1/N 2 expansion called topological expansion (N is the size 
of the matrix): 

CO 

F = Y^N~ 2h F {h) (1.1) 

h=0 

The authors of invented an efficient method to compute recursively all the F^'s, 
and they improved it in 0. 

The 2-matrix model is conjectured to have a similar 1/N 2 topological expansion. 
Indeed, this expansion is the main motivation for applications to 2-dimensional quan- 
tum gravity jTj^, |25|, because each is the partition function of a statistical physics 
model on a genus h surface. There is at the present time no rigorous prooiQ of the 
existence of such an expansion, but, assuming that it exists, the aim of the present 
work is to give a method to compute recursively the terms of the expansion, similar to 
that of II. 



The 2-matrix model was first introduced as a model for two-dimensional gravity, 
with matter, and in particular with an Ising field P3LO]. The diagrammatic expansion 



2 The Riemann-Hilbert approach seems to be the best way to prove the existence of the 1/N 2 
topological extension as in |22]. The Riemann-Hilbert problem for the 2-matrix model has been 



formulated |10|,|9[ pi B2| , and seems to be on the verge of being solved 
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of the 2-matrix model's partition function is known to generate 2- dimensional statistical 
physics models on a random discrete surface [p^,p!7|,p3|: 

N 2 F = -\nZ = e- Action (1.2) 

surfaces matter 

where the Action is the matter action (like Ising's nearest neighboor spin coupling) plus 
the gravity action (total curvature and cosmological constant) |19|| . The cosmological 
constant couples to the area of the surface, and N (the size of the matrix) couples to 
the total curvature, i.e. the genus of the surface. The large N expansion thus generates 
a genus expansion: 

CO 

F = ^N' 2h F w (1.3) 

where F^ h ' is the partitrion function of the statistical physics model on a random 
surface of fixed genus h. 

p(h) _ e -Action ^ 

genus h surfaces matter 

The leading term F^ is the planar contribution. It can be computed by many different 
methods, for instance the saddle point method (along the method invented by |35j and 
rigorously established by |3~0|D or the loop equation method which we will explain 
below. Our goal in this article is to compute F^ and present an algorithmic method 
for computing F^ h ' for h > 1. We generalize the method of ||. 

1.1 Reminder: the one- matrix model 

For a given polynomial V{x) of degree d + 1, called the potential, and a given integer 
N, we define the partition function Z and the free energy F as: 

Z = Q -N*F = J dMe -NtrV(M) ^ 

where the integral is over the set of hermitian matrices of size N, with the measure 
dM equal to the product of Lebesgue measures of all real components of M. 

The free energy has a 1/N 2 power series expansion (under the 1-cut asumption, 
this will be discussed in more details in section |4.1|) : 



F = Y^N~ 2h F {h) (1.6) 

h=0 

In their pioneering work, the authors of invented a method to compute recur- 
sivelly all the F^'s. In particular they found that the genus one free energy for the 
1-matrix model is: 

F« = -l In ((6 -a) 4 M(a)M(6)) (1.7) 
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where b, a and the polynomial M(x) are related to the large N limit of the average 
density of eigenvalues of the matrix: 

p(x) = ^-M(x) v / {x - a) (b - x) , M(x) = Pol V'{x)/y/{x - a)(x - b) 

(1.8) 

This density has a compact support [a, b], and M(x) is a polynomial of degree d — 1 
(the 1- matrix case is re-explained in section |7|) . 

We are going to extend this kind of expression (eq. ( |1 . 7| ) ) for the 2- matrix model. 

1.2 Outline cof the article 

• In section ^ we introduce the definitions and notations, in particular we define 
the 1-loop functions and 2-loop functions, and the loop-insertion operators. 

• In section ||] we explain the loop equation method, and derive the "Master loop 
equation" . 

• In section U we observe that, to leading order, the master loop equation is an 
algebraic equation of genus zero, and we study the geometry and the sheet- 
structure of the underlying algebraic curve. 

• In section |5] we use the algebraic master loop equation to compute all the previ- 
ously defined loop functions to leading order. 

• In section ^ we include the previously neglected 1/N 2 term in the loop equation, 
and we compute the 1-loop function Y(x) to next to leading order. Then we 
derive the next to leading order free energy by integrating Y^(x). We also 
discuss how to compute higher order terms. 

• In section [7| we compare with the results previously known for the one matrix 
case. 

• section E is the conclusion. 



2 The 2-matrix model 
2.1 The partition function 

For two given polynomials V\(x) of degree d\ + 1 and V 2 (y) of degree d 2 + 1, called 
the potentials, and a given integer N, we define the partition function Z and the free 
energy F as: 

Z = e~ N2p = j dMidM 2 e~ NtI [Vi(M 1 )+v 2 (M 2 )-M 1 M 2 ] ^ 
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where the integral is over the set of pairs of hermitian matrices Mi and M 2 of size 
N, with the measure dMidM 2 equal to the product of Lebesgue measures of all real 
components of Mi and M 2 (see for definitions). 
The potentials V\ and V2 are written: 

d\ rig * 



-X 



W = E 



fc + i ' v ^ ^k + i 

k=0 k=0 



y 



fe+i 



(2.2) 



Remark: eq. (J2T 
i.e. when: 



is well defined only when the potentials are bounded from below, 



[d x and d 2 are odd g dl+1 > 0, g* d2+l > 0} (2.3) 

However, it is known (l9j how to give a meaning to |2.1| for potentials non bounded from 
below, in the large N limit (somehow the tunnel effect allowing the matrices to "escape 
from the potential wells" is suppressed like 0(e~ N ) in the large N limit, it plays no role 
in the 1/N 2 expansion). Therefore, we will consider arbitrary polynomial potentials 
which do not necessarily satisfy condition (|2.3|). For simplicity, we will assume that 
the potentials V\ and V 2 are realQ. 

In addition, we will assume that the potentials V\ and V2 are generic, i.e. they are 
not critical (we will explain what is a critical potential in section |4.4.1| , see ||16|| ). 

The partition function F can be interpreted (through the Feynmann graphs expan- 
sion, see []B5|,[25|]) as the generating function of surfaces made of polygons with a spin: 
spin "up" (+) for the polygons generated by the vertices of M 1; "down" (-) for the 
polygons generated by the vertices of M 2 . 



where £ is the ensemble of surfaces made of up-polygons with at most d\ + 1 sides and 
down-polygons with at most rf 2 + 1 sides. Aut(S') is the group of automorphisms of the 
graph of edges of S, h(S) is the genus of S, n k (S) (1 < k < d\ + 1) is the number of 
up fc-gons of S and fij(S) (1 <j < d 2 + 1) is the number of down j-gons of 5*. 
We write digrammatically: 




3 The derivation of the loop equations in section || should be slightly modified in the case of non-real 
potentials, because some change of variables we consider would be no longer hermitian. But most of 
the results derived in this article are still valid for complex potentials. 
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This is why the two-matrix model can represent the Ising model on a random 
discretized surface [^,|T^], and other 2d statistical physics models on a random surface 
carrying some type of matter It is clear from eq. ( |2.4j) that relevant applications 
to 2d statistical physics on a random surface involve potentials which violate condition 
(all the g^'s and gi§ must be negative in order to have positive weights). 



2.2 Definition of the loop functions 

We introduce the following expectation values computed with the probability weight 

i p -JVtr [V 1 {Wh)+V 2 {M2)-M 1 M 2 \ . 



1 

N 



tr M^M l 2 ) 



(2.5) 



We will be particularly interested in their large N limits and large N expansion. 

Tkj can be interpreted as the partition function of a statistical physics model on 
surfaces with one boundary of lenght k + 1 (k contiguous + and j contiguous -): 



Tk,, 




+ 



It is more convenient to introduce generating functions which contain all the T^'s 
at once: 



2.2.1 One-loop functions 

We define the following formal functions of the complex variables x and y: 



W(x) = — ( tr — -— — 



W(y) = — I tr 



k=0 



k+l 



oo rp 

yk + 1 

k=0 y 



Y(x) = V((x) - W(x) (2.6) 



X(y) = V&y) - W(y) (2.7) 



These functions are formally defined only through their large x and large y expansions, 
they are only a convenient rewritting of the collection of T^q and To,&- The sum is not 
necessarily convergent. 

We will find below, that W^(x) = lim A r^ 00 W / (x) (resp. W(y)) satisfies an algebraic 
equation, and thus W^°\x) (resp. W^°\y)) is an analytical function of x (resp. y), in 
the complex plane, with a cut. The radius of convergence of eq. ( p.6| ) is finite in the 



large N limit. Their analytical structure will be dicussed in details in part 4.4 
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The functions W(x) and W(y) are called the resolvents, they provide most of the 
information about the statistical properties of the spectra of Mi and M 2 . For instance 
the location of the cut of W^\x), is the support of the large N average density of 
eigenvalues of Mi, and in the large N limit, the average density of eigenvalues is 
proportional to the discontinuity of W^\x) along its cuts. 

Their diagrammatic expansion is the generating function for statistical physics mod- 
els coupled to gravity, on a surface with one boundary (a disc if genus zero), this is 
why they are called one-loop functions. 



w( X ) = j2*- h - 1 E E 



— Action 




+ 



(2- 



surfaces, matter 

boundary — l-± up 



where the sum carries on all discretized surfaces with one boundary made of li spins 
+. 



—Action 




+ ... (2.9) 



«2=0 



surfaces, matter 

boundary — Z2 down 



where the sum carries on all discretized surfaces with one boundary made of l 2 spin -. 
We also define: 



W(x, y) — -J- I tr — — : — — 

v ' yj N \ x-M iy -M 2 



EE 



T, 



k.i 



(2.10) 

k=0 1=0 

This function provides information on the correlation between the spectrum of Mi 
and that of M 2 . Its diagrammatic expansion generates surfaces with one "bi-colored" 
boundary, i.e. a boundary of lenght li + l 2 made of li consecutive spins up, followed 
by Z 2 down: 



w(x, y ) = j2J2 x ~ h ~ l v~''~ 1 



E E 



-Action 




l 1= 0l 2 =0 



We also define: 



surfaces, matter 

bnd = l\ up + I2 down 



N \ x- Mi y - M 2 



(2.11) 
(2.12) 
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u (x , y) = LUYM^W. i \ (2 . 13) 

iV \ x - Mi y-M 2 / 

which are polynomials in one of the variable and analytical near oo for the other 
variable. 

And we define 

N \ x - Mi y - M 2 / V ; 

which is a polynomial in both variables, of degree (c?i — 1, d 2 — 1). 

In appendix C we compute another one-loop function H(x, y, x f , y'), which generates 
surfaces with one quadri-colored boundary, i.e. surfaces whose boundary is +-+-. 
Somehow it is the random surface analogus of a square with opposite sides of the 
same spin and adjacent sides of opposite spin. It maybe very usefull to compute that 
kind of function to study how spin percolates from one side to the other. It may 
have applications to compute boundary operators in 2d statistical physics on a random 
surface. 

2.2.2 Two-loop functions 

We also need to define the following correlation functions, which contain the product 
of two traces. Their diagrammatic expansion is the generating function for statistical 
physics models coupled to gravity, on a surface with two boundaries (cylinder-like), 
this is why they are called two-loop functions. 



n ^^=S = ( te ^ te ^ > - jvaw(acWa ° ' (2 - l5) 



dU(x,y) / 1 V 2 \y)-V>(M 2 )^ 1 x - 



U ^ = -W0f = \ tr ,-Mi y-M 2 " tr— , 

(2.17) 



dU(x,y) I V[{x)-V[{M X ) 1 _ 1 x - 



(2.18) 
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P(x,y;x') 



dP(x,y) 
dVi{x') 



tr V({x) - V{(Mi) Vi{y) - V 2 \M 2 ) ^ 



1 



x - Mi 
-N 2 P(x,y)W(x') 



y-M 2 



x' - Mi 



where the "loop- insertion" operators (see 0) are formally^ defined as: 







dVAx') 



k d 
\ x' k dg k 







dV 2 (y>) 



k d 



(2.19) 



[2.20) 



their action on the partition function or on expectation values inserts a trace, i.e. a 
boundary (a "loop") for the ensemble of surfaces. For instance we have: 



W(x) 











dVAx) 











oVi[x') oVi{x) oVi{x') 

(2.21) 

The diagrammatic expansion of F gives surfaces with no boundary, W(x) gives surfaces 
with one boundary, Q(x,x') gives surfaces with two (spin up) boundaries, and Ct(x, y) 
gives surfaces with one spin up boundary and one spin down boundary. And so on: 
inserts a + loop, inserts a - loop. 



Q(x; x') 




Q(y;x') 




3 The loop equations 
3.1 Generalities 

The loop equations are a consequence of the reparametrization invariance of the matrix 
integral. Basically, one writes that under an infinitesimal change of variable of the type: 

M — .\/| M • e/(Mi, M 2 ) (3.1) 

the partition function is left unchanged to order 1 in e. We have to compute the 
variation of the measure and the variation of the potential: 

Z = J dMidM 2 e~ Wtr[yi(A ^ l)+F2(M2) ^ lM21 

4 By formally, we mean they make sense only through their large x (resp. y) expansion. If k is 
larger than d±, d/dg^ means: take a potential of degree > k, and then take the derivative w.r.t. <?j, at 
g.j = for j >di. 
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dM 1 dM 2 (1 + eJ(M x , M 2 )) e~ Ntt m^+v^M^-MyM^]^ _ e ^( Ml , M 2 )) 
(3.2) 

where J (Mi, M 2 ) is the jacobian of the change of variable computed to order 1 in e: 

det = 1 + eJ ( M i> M 2) + °( e2 ) ( 3 - 3 ) 

and K(Mi, M 2 ) is the variation of the action to order 1 in e. We must have: 

(J(Mi,M 2 )) = (K(Mi,M 2 )) (3.4) 

We are going to give the receipeQ to compute J (Mi, M 2 ) for various change of variables 
f(Mi,M 2 ). 

Remarks: 

- It is important to take into account the non-commutativity of matrices. 

- It is also important to choose a change of variable which is a hermitian matrix. 

- We are presenting the method for variations of Mi but of course the same applies 
to M 2 and we will use both changes of variables in the body of this article. 

3.2 split and merge rules 

We will be interested in change of variables of the following formQ: 

• f{M li M 2 )=A^ r B: 

where the matrices A and B are functions of Mi and M 2 . The variation of the 
measure is: 

J(Mi,M 2 ) = tr [A - ) tr f - B ) + contributions from A, andB. 

\ x — Mi J \x — Mi ) 

(3.5) 

somehow, you split the trace into two traces whenever you meet a X ^ M term 
outside a trace, and each side receives a — \r- 

' x— Mi 

• f(Mi,M 2 )=Atr (B^): 

where A, B are functions of M\ and M 2 . The variation of the measure is: 

J(Mi,M 2 ) = tr ( A B 1 + contributions from A, andE. (3.6) 

v ; V x-Mi x-MiJ ( J 

somehow, you merge the traces whenever you meet a x \ Ii term inside a trace, 
and duplicate the X 2 M term (one before B, one after B). 



5 Easy to prove directly order by order in x. We leave the proof to the reader. A similar method 
is used in |2(|. 

-Mi 



6 As before, X _^ M = J2k x k 1 M\ is a "formal" notation for the collection of My , k = . . . oo. 
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• The variations of A and B are computed recursively by the chain rule. 

The variation of the action is simply: 

^K(M U M 2 ) = tr {Vi{M 1 )f{M ll M 2 ) - M 2 f(M u M 2 )) (3.7) 
and we will use the following trick very often: 

V 'x-M x ) 1( J \x-Mx J V a; -Mi J y ' 

where it is important to notice that the second term is polynomial in x. 

3.3 Factorization "theorem" 

Then, we will make the asumption that in the large N limit, the average of product of 
traces decouple: 

i ( tr A tr B) ~ i ( tr A) i ( tr B) + 0(1/N 2 ) (3.9) 
more precisely that the connected correlation function 

(tr AtrB) c = (tiAtrB) - (tr A) (tr B) (3.10) 

has a large N limit. 

The factorization "theorem" is not a theorem, it was never proven for the 2-matrix- 
modelQ, but there are plenty of heuristic reasons to believe that it holds, and it was 
used by nearly all authors in this field. As we said in the introduction, we are not going 
to prove the existence of the 1/N 2 expansion, but assuming it, we are going to give an 
algorithmic method to compute it. 

3.4 The master loop equation 

Consider the change of variable: 

M -i — M i ■ < ' 



x - Mi 



which is indeed a hermitian matrix. It leads to (there is no variation of the measure 
since 5M 2 is independent of M 2 , the RHS is the variation of the potential): 

= ( tr l —VXM 2 ) (3.11) 



7 It was proven rigorously for the one-matrix model in 
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i.e. 



tr 



V 2 (M 2 ) ) = xW(x) - 1 



N \ x - Mi 

Now, consider the change of variables (which is indeed a hermitian matrix): 



(3.12) 



.Ui • - M ■ < 



V±{y) - V 2 \M 2 ) 1 



1 V 2 \y) - V 2 \M 2 ) 



y-M 2 x-Mi x-Mi y-M 2 

it leads to (the LHS is the variation of the measure, the RHS of the potential, and we 
use the cyclicity of the trace): 

1 1 V>{y)-V 2 \M 2 ) 



tr 



tr 



x - Mi x - Mi y - M 2 
.V { tr Vl{Ml ) V * iy) ~ ^ {M2) tr 1 



x - Mi y - M 2 
i.e. (using eq. fl3.12Q ): 



x - Mi 



-M 



y-M 2 



{y - Y(x))U(x, y) = V 2 \y)W{x) - P(x, y) - xW(x) + 1 - — U(x, y- x) 



(3.13) 
(3.14) 

(3.15) 



Notice that U(x,y) is a polynomial in y, therefore it is finite for y = Y(x). If we 
choose y = Y(x) we get (using W(x) = V{(x) — Y(x))\ 



V = Y(x) 

(V 2 \y) - x)(V{(x) -y)- P(x, y) + l = ^U(x, y- x) 



(3.16) 



which we write: 



E(x,Y(x)) = —U(x,Y(x);x) 



(3.17) 



where E(x,y) is a polynomial in x (degree d\ + 1) and y (degree d 2 + 1): 



E{x,y) = {V{{x) - y){V 2 \y) - x) - P{x,y) + 1 



(3.18) 



eq. ( 3.1 7|) is called the "master loop equation" [41], we will see below that it allows to 
determine the function Y(x) and all other functions in the problem. 

3.5 1/N 2 expansion 

We assume that the function Y(x) can be written as a power series in 1/N 2 : 



Y(x)=Y(°\x) + -^-Y( 1 \x) + ... 



(3.19) 



Such an expansion is believed to hold in the "1-cut case", which will be explained below 
(section |4"TD . 
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To leading order, the loop equation is an algebraic equation: 

£ (0) (x,y (0) (a;)) = (3.20) 

where 

E^(x,y) = {V{{x) - y){V 2 \y) - x) - P {0 \x,y) + 1 (3.21) 

The polynomial P(°\x, y) has degrees (d\ — 1, d 2 — 1), it is determined by the large 
x and large y behaviors of Y (x) and X(y), as well as some analytical requirement, 
namely the "one-cut asumption". As we will see below, in the 2-matrix model, 
this asumption is replaced by the asumption that the algebraic curve E(x, y) = 
has genus zero. In the next section we will study this algebraic curve, and we 
will explain how to compute the polynomial P^°'(x,y). 

Let us note that this algebraic equation was derived by amny authors and various 
methods in the literature: withe the loop equation method, the existence of an 



algebraic equation is proven (and computed for small degree potentials) in [f|l 
, and was presented for arbitrary polynomial potentials in the appendix of 



and P^| . The authors of |45j and |R| derived it by the "saddle point method" 



inspired from [p^| , and from the 2-toda hierarchy dispersionless limit, and that 
algebraic equation was used by many other authors in particular cases. And 
notice that it is the same equation as the spectral curve of any of the 4 dual 
differential systems in ||. 

To order 1/N 2 we expand: 

E{x, y) = E®(x, y) - ~^P {1) (x, y) + . . . (3.22) 

where P^(x,y) has degrees (d\ — l,d 2 — 1). We will see how to compute it in 
section 0. 



We easily get: 



Y 0) M _ U {0) (x,Y(°\x);x) + PW(x,YM(x)) 

[X) ~ d y E(°)(x,Y(°)(x)) {6 - Z6) 



In this expression, all terms in the RHS are computed to leading order, except 
the polynomial P^(x,y). But P^\x,y) is completely determined (in terms of 
leading order functions only) by the condition that all the unwanted zeroes of 
the denominator (those which don't correspond branch-points) cancel. In other 
words, once we know the functions Y(x) and U(x,y;x') to leading order, we 
can compute Y(x) to order 1. Similarly, we can compute every loop function to 
order 1. 
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• We can repeat the procedure to find higher orders. Once we know all the loop 
functions up to order h, we can compute them to order h + 1. Of course the 
expressions become more and more complicated with increasing orders, but the 
method is the same. We will illustrate our method by computing the order 1 
only. 

• once we know Y^ h '(x), we can find the order h free energy F^ h \ by integration: 

= -dF^/dVt(x). 

4 Leading order: algebraic equation 

While we study the leading order, we will drop the subscripts (0) for lisibility. 

4.1 Reminder: one-cut asumption for the 1-matrix model 

In the 1-matrix model, the existence of the 1/N 2 topological expansion is related to 
the 1-cut asumption. We refer the reader to [|1^,[18],|22| for more detailed explanations, 
and we just summarize the main idea. 

When there is more than one cut, i.e. when the density of eigenvalues has a dis- 
connected support, there are tunnel effects between the different cuts which generate 
an oscillating function of iV (a theta- function). Such a periodic function cannot have 
a 1/N 2 expansion. The frequencies of oscillations are the filling fraction of eigenvalues 
in the cuts, in other words, the frequencies are equal to the integrals of W(x) along 
the non-trivial cycles of the hyper-elliptical surface generated by (x, W(x)) (In the 1- 
matrix case the surface is hyperelliptical, i.e. W(x) obeys a quadratic equation). It is 
only in the one-cut case that there is no oscillations, and the 1/N 2 expansion can exist. 

We expect the same idea to hold for the 2-matrix model. The method of JOS can 
be adapted to the 2-matrix case: the frequencies are the contour integrals of W(x)dx, 
which are the same as the contour integrals of Y(x)dx (which are also the same as 
the integrals of X(y)dy by integration by parts), along the non-trivial cycles of the 
algebraic curve (x,Y(x)). The 1/N 2 expansion cannot exist if there are oscillations. 
Therefore we will require that there is no oscillation, i.e. that there is no non-trivial 
cycle, i.e. that the algebraic curve has genus zero. 

Notice a difference between the 1-matrix case and the 2-matrix case: 

The number of connected parts of the support of the density is actually the number 

of cuts in the first sheet only. There may be cuts in other sheets. 

- In the 1-matrix case, the curve is hyper-elliptical, it has only two sheets, thus 

the first sheet contains all the cuts. The number of non-trivial cycles is equal to the 
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number of cuts minus one. Therefore in the 1-matrix model, the genus zero asumption 
is equivalent to the 1-cut asumption and is equivalent to a density with connected 
support. 

- In the 2-matrix case, the curve is only algebraic. It can have more than two sheets, 
so that the non-trivial cycles are not necessarily in the first sheet. One cut in the first 
sheet is not equivalent to genus zero. In fact, there are g?2 + 1 sheets and at least d<i cuts, 
g?2 — 1 of them are NOT in the first sheet. In the 2-matrix case, the connectedness of 
the support of the density is a necessary condition, but is not sufficient. The necessary 
and sufficient condition is the genus zero condition. 

Therefore, in the 2-matrix model, we claim that the existence of the 1/N 2 expansion 
is related to the asumption that the algebraic curve E^(x,y) = has genus zero. 

From now on, we will consider only the genus zero case. We would like to emphasize 
that many of the results derived below for the leading order loop functions, can be 
derived for arbitrary genus with no difficulty (for instance the functions Q are always 
the Bergmann kernels), but are out of the scope of this article. Many of them are 
presented in |J. 



The previous discussion applies to the case where the partition function eq. Q2.1| ) is 
well defined, i.e. the potentials satisfy condition eq. ( fOD . However, there is another 
case of "existence" of the 1/N 2 expansion: when the genus h free energy is formally 
defined by its diagrammatic expansion eq. ( |1.4j ) and eq. ( |2.4| ). The formal series for 
F^ h ' is obtained by expanding in the vicinity of a gaussian integral [19, 14], i.e. by 



considering that the potentials are "small deformations" of quadratic potentials. When 
both potentials are quadratic, the algebraic curve E(x, y) = clearly has genus zero. 
The formal expansion in power series of the g^s and #*'s can not change the genus. 
Therefore, the formal diagrammatic expansion of F^ h > is also obtained by a genus zero 
asumption. 

4.2 Genus zero asumption 

Since we require that E(x,y) = be a genus zero algebraic curve, there must exist 
a rational uniformization. Such a rational uniformization appears naturally in the 



framework of bi-orthogonal polynomials as shown in p, £3], [24], [U| . That leads us to 
choose it of the following form: 

x = X{s)= 1 s + Y j ^ , y = y( 8 ) = l + J2p ja i (4.1) 

k=0 j=0 
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This means that for every (x, y) which satisfy E(x,y) = 0, there exists at least one 
s such that x = X(s) and y = y(s). 

We note the inverse functions: 

x = X(s) <-> s = a{x) , y = y(s) <->• s = a{y) (4.2) 

The functions a(x) and a(y) are multivalued, we will discuss their sheet structure 
below. The functions Y(x) and X(y) are: 

Y(x) = y(a(x)) , X(y) = X(a(y)) (4.3) 

They are multivalued too, and their sheet structure will be discussed below. 

4.3 The parameters oik and j3j 

The a's and /3's are merely a reparametrization of the coefficients of the potentials. 

Indeed, since for large x we have: V{(x) — Y(x) = W(x) ~ - + 0(1/ x 2 ) (this 
corresponds to s — > oo), and for large j/ we have: V 2 '(?/) — X(y) = W(y) ~ - + 0{l/y 2 ) 
(s — > 0), we must have: 

W - ~™ -- + 0(l/ s 2 ) 

~ s ^o -~ + 0(s 2 ) (4.4) 

7 

Both these two equations imply that (compute Resy(s)A"(s)ds = — ResAf(s)3^'(s)ds): 



min(di ,^2) 

7 2 = -1 + ^ ka kPk 
fc=l 



(4.5) 



And eqs.(fOj) allow to write the coupling constants gi and g\ as functions of the 
ctfc's and /Vs (compute 'ResV{(x)x l ~ 1 dx and ResV 2 '(?/)?/~ :l d?/): 



dx-l-l I d 2 [rfi/(p+l)]\ n , v ; </.. / <t x ,,. 



j dl-l-l / d 2 [Oi/y>-t-J.jj \ (; + V i )! 

^+i=-3k s n s * /! n i f a ^(i+k+i+E P (p+^p) n v - 

' fc=-i \p=o i P =o / '■llpV p =o V 7 

(4.6) 

^ (fe-i-l / d a [<fe/(p+l)]\ ^ + i p y d i / a \h> 

e n s * /! n P z f ^ a c+^+i+E P ( P +i)v) n ( ) 

7 fc =_i y P =o i p =o / '-iip'p- P =o v 7/ 

(4.7) 

where we take the convention that (3j — if j > d\ and atj = if j > d 2 , and 
a_i = = 7- 
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As an example, for / = d\ and I — d\ — 1, the formula reduces to: 

9*+i = ^k ' ^ = ^zr - diao^ , ••• (4-8) 

In principle, it should be possible to revert these formula, and compute the a's and 
/3's as functions of the coupling constants. This can be done at least numerically. 
Remark: 

Notice that there might exist more than one solution. Only the solution which leads 
to an absolute minimum of the leading order free energy should be considered. 
If the potentials V\ and V 2 are real (which is the case for most physical models), the 
density of eigenvalues of the first (resp. second) matrix is the imaginary part of Y(x) 
(resp. X(y)) along the cut, and the density must be positive. Any solution for the a's 
and /3's which does not satisfy that is unacceptable. If there is no acceptable solution, 
this means that our choice of V\ and V 2 does not correspond to a genus zero case. In 
general, if V\ and V 2 have only one well, or if they are close to quadratic potentials, we 
have a genus zero solution. 

4.4 Sheet structure of the function Y(x) 
The function a(x) is multivalued, indeed the equation 

X(s) = x (4.9) 

has d 2 + 1 solutions, which we note 

{(T (a;),(7i(a;),...,(T d2 (x)} (4.10) 

and therefore, the function Y(x) = y(a(x)) is multivalued with d 2 + 1 values which we 
note: 

{y (z),n(z),-..,^ 2 (z)} where Y k (x) = y(a k (x)) (4.11) 

We define cr(x) = <j (x) to be the "physical sheet" solution, i.e. the one such that 
cr(x) — > oo when x — * oo, i.e. the one for which W(x) = V[(x) — Y{x) ~ 0(1/ x). The 
other solutions a k (x) with k — 1, . . . , d 2 all tend to zero when x — > oo. 

4.4.1 Endpoints and cuts 

The functions Y(x) and a(x) have cuts, whith a generic square root behaviour near the 
endpoints, which means that the derivative with respect to x diverges at the endpoints. 
since Y'(x) = p^, the endpoints are solutions of s = or solutions of the equation: 

X'(s)=0 (4.12) 
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which has di + 1 solutions. We note them: 



{e ,ei,...,e d2 } (4.13) 

the endpoints are thus located at the X(ek)'s, as well as at infinity. In the physical 
sheet, the function Y(x) is analytical at infinity, and the only cut in the physical sheet 
is a line joining two endpoints of the type X{e^). In the other sheets, the function 
Y(x) behaves as x 1 ^ 2 near infinity, the cuts are lines joining some X(ek) to oo. There 
are exactly d 2 cuts, in agreement with the fact that the genus is zero. 

Similarly, the function X(y) has exactly d\ cuts, whose endpoints are y{ej)'s and 
oo, where 

{e ,gi,...,g dl } (4.14) 

are the di + 1 roots of 

y'(s) = . (4.15) 

The cut in the physical sheet is a line joining two endpoints of the type y(ek), and the 
cuts in the other sheets are lines joining some y(e~k) to oo. 

As said in the introduction, we assume that the potentials are not critical, which 
means that the algebraic curve has no singular points (other than s = and s = oo). 
That means that all the and e\ must be distinct (examples of critical potentials are 
given in |TB[]: those which produce the rational minimal conformal models (p,q))- The 



non-criticality asumption, is crucial for the derivations we are going to present below: 
we will often have X"[e^) or y\e^) in the denominators. 

4.4.2 Sheets in the s-plane 

We define the x-sheets (resp. y-sheets) in the s-plane, as domains which contain only 
one crfc(x) (resp. <7k(y)) when x (resp. y) sweeps the complex plane. In other words, 
the function X(s) (resp. y{s)) restricted to one domain is one to one. 

The domains are separated by closed curves in the s-plane, which are the projections 
of the cuts in the x— plane (resp. y-plane). There is some arbitrariness in the choice of 
domains and cuts (one could take the cuts as straight lines or any other curves which 
link the endpoints). The only constraint is that all the cuts go through the endpoints. 

The x-sheets (resp. y-sheets) have the following properties: 

- The physical sheet is the sheet which contains oo (resp. 0). 

- All the d 2 (resp. di) other sheets contain (resp. oo). 

And the cuts: 

- The cut in the physical sheet corresponds to a contour in the s-plane which encircles 
0. If the potentials are real, we expect the cut in the x-plane (reps, y-plane) to be on 
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the real axis, and its s-plane image to be a contour which crosses the real axis, and 
symmetric with respect to the real axis. The two endpoints are solutions of X'(s) = 
(resp. y'{s) = 0). 

- The other d 2 — 1 (reps. d± — 1) cuts are contours in the other sheets. They all meet 
in (resp. oo). One of their endpoint is a solution of X'(s) = (resp. y(s) = 0) and 
the other endpoint is at s = (reps, s = oo). 

notational remark: when we write a(x) this means that a(x) is considered as a 
multivalued function. For instance the multivalued function Y(x) is equal to 

Y(x)=y(a(x)) (4.16) 

We write <Jk{x) or Yk(x) only when we want to specify in which sheet is x. By default, 
we will consider that x is in the physical sheet, and therefore <7q(x) is in the domain 
which contains s = oo, and <7k{ x ) with k > 1 is in a domain which contains s = 0. 
Note that any symmetric function of (cr (x), . . . , (Jd 2 (x)) is a regular function of x, i.e. 
it has no cut. 



4.5 The polynomial P(x, y) 

The polynomial P(x,y) is found from E(X(s),y{s)) = 0, i.e. : 

p(x( s ),y(s)) = (v;(x(s)) - y(s))M(y(s)) - x( 8 )) + 1 (4.17) 

both sides are Laurent polynomials of s of the same degree. By indentifying the 
coefficients of powers of s on both sides, one can compute explicitely all the coefficients 
of the polynomial P(x,y). The expression of P(x,y) in terms of and f3j can be 
written explicitely with the determinant formula: 

/ 7 a — x cci ... aa 2 ... \ 

••. ■•. 



E(x,y) 



di+d.2 



det 






V 



7 ao — x Oil 

(h Po-y 7 



«d 2 








Pi Po-y 7 / 



(4.18) 

This determinant of size d\ + d 2 + 2 is the discriminant which vanishes if and only if 
the Laurent polynomials in s: X(s) — x and y(s) — y have a common root. 
There is also the formula: 

E(x,y) = -(g dl+ i) d2+1 f[f[(l<Ti(x) ~ itM) (4-19) 

—0 j=0 
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indeed, it is a symmetric function of the 0^(2) and &j(y), with no poles, therefore it is 
a polynomial in x and y, and it vanishes when X(s) — x and y(s) — y have a common 
root. 

Remark: Had we not made the genus zero asumption, we could determine P(x,y) by 
a determinant formula^ very similar to eq. (|4.18 ). We can also determine P(x,y) by 



requiring that E(x,y) = is a genus g curve, and that all the b-cycle integrals of 
Y(x)dx vanish. There is one (or more) solution for all g. For any g between and 
|" rii<fe-i j we can compute all the loop functions (to leading order), and therefore we can 
compute the free energy F^°\g) as a function of the genus g. The actual value of g is 
the one which minimizes F^°\ The genus zero asumption is thus that the minimum is 
obtained for g = 0. 

5 Leading order loop functions 

5.1 Some leading order Loop equations 

The change of variables (we drop the 0(1/N 2 ) terms): 

• Mi -»• Mi + + implies: 

(y - Y(x))W(x, y) = W(x) - U(x, y) = V{{x) - Y(x) - U(x, y) (5.1) 

• M 2 - M 2 + + implies: 

(x - X(y))W(x, y) = W(y) - U(x, y) = V^y) - X(y) - U(x, y) (5.2) 

. M, - M, • 6 \ f V ^ y) -f Ah) + e «zffiH) 1 implies: 
1 1 x—Mi y—M2 y—M2 x—M\ 1 

(y - Y(x))U(x, y) = -P(x, y) + V^y)W(x) - -J- < tr -A— K 2 '(M 2 ) > (5.3) 

N x — Mi 

. M 2 ^ M 2 + e MEhgMl 1 + e^ SMzgW imp i ies: 

z z 1-M1 y—M-i y—M2 x—M\ r 

{x - X{y))U(x, y) = -P(x, y) + ^(x)^(y) - 1 < tr — L^Mi) > (5.4) 



M 2 -> M 2 + impl 



les: 



4 < tr —-VUM*) >= xW(x) - 1 (5.5) 

N x-Mi 



Mi -> M x + implies: 



8 Thanks to J. Hurtubisc for that remark. 
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5.2 The function U(x,y) 



Using eq. 
simply: 



and eq. Q5.5|), we find that to leading order, the function U(x,y) is 



U(x,y) = V^y)-x + 



E(x,y) 
y ~ Y(x) 



(5.7) 



Notice that this function has no pole, and it is indeed a polynomial in y. 



5.3 W(x,y) 

Using eq. (|5~2] ) and eq. Q5.71 ), we find that, to leading order, the function W(x,y) is: 



(5i 



W(x,y) = l- 



E(x,y) 



[x-X(y))(y-Y(x)) 



It can also be written: 

W(x,y) = l 



? dl Y 2 <7o(*)*o(v)* f[f[ {(Ti{x) 



7 ) a Q (x) -a (y) fj-f^ 



(5.9) 



it has a pole when a^x) = cr^y). 

We are now going to compute the two-loop functions. 

5.4 The functions Q(x;x') and Q(y;x') 



These functions were first computed by |UJ using the orthogonal polynomial's technics, 
they can also be found by another method, namely: they must be rational functions 
of a(x) (resp. &(y)) and cr(x'), and they are determined completely by their poles and 
behaviours near oo. The result is the following: 

1 



Q(x; x') 



dVi (ar 



-W(x) = d x d x > In ■ 



+ 



x — x' 
1 



+ 



a'(x)a'(x') 



(x-x') 2 X'(a(x))X'(a(x'))(a(x) - a(x')) 2 
and since Y(x) = V((x) — W(x) we have: 



(x — x') 2 (cr(x) — a(x')) 2 

(5.10) 







dV x (x') 



-Y(x) 



a'(x)a'(x') 
(cr(x) - a(x')y 



-d x d x > In (a(x) - a(x')) 



(5.11) 



i.e. » d Y(x) is the Bergmann kernelf]. 



dVi(x') 

9 Without the genus zero asumption, one has dv ^ x ,^ (x) — Bergman kernel, but dyf^') ^( a 
square of the Szego kernel +0(1/N). 
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At x — x' we get the Schwarzian derivative of cr{x) 



n(x , x)= ^.\qfi ( , 12) 

6 a'(x) 4 a'(x) z 



Similarly the function Q(y;x') is: 



fi <* = wvf^ = 111 cm - = -j£0§fr- (5 ' 13) 

am™ = = d ' d " 1,1 - a{x ' ]) (5 - 14) 

5 More loop equations 

• The change of variables 5M 2 = x } Ml tr X ,^ M implies: 

{*Thr*^*hr)r' a{xi ' r) (5 ' 15) 

• Using this result, the change of variables 5M\ = — 1_ Y^Mh^Mlt i x 1 -f h.c. 
implies: 

(y-Y(x))U(x,y;x') = (V^y) - x - U(x, y))Q(x; x') - P(x, y; x') 

_ dx u(^-uj*',y) (5 . 16) 

i.e., using eq. (|5.7|): 

(y -Y(x))U(x,y;x) = -tt(x;x) — - -P(x,y;x) - d x , 

y — Y (x) x — t 

(5.17) 

which can also be written, using eq. (|5.7|) again: 



(y-Y(x))U(x,y;x') = -(n(s; x>) + J^y)^^ ~ ^(wO 

x (x — x')(y — Y(x'j) 



The change of variables 8M\ = y \ l2 tr x ,* Mi implies: 

1 

' y-M 2 ~~ x'-Mt 



XxV^Mx)—^— tr — L- ■ ) = 2 /fi( 2 /;x')-^^(x , ,i/) (5.19) 
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• Using this result, the change of variables 5M 2 = X 1 }^ — Vr tr tV + h.c. 
implies: 

(x - X(y))U(x, y; x') = Q(y; x'){V({x) - U(x, y) - y) - P(x, y; x') + d x ,W{x', y) 

(5.20) 

which can also be written: 

(x - X(y))U(x, y; x') = -P(x, y; x') - Q(y; x') E ^' y \ + d x ,W(x', y) (5.21) 

x - X(y) 

5.6 Determination of the Polynomial P(x,y;x') 

In particular, if for any k = 0, ...,c?2, we choose y = Y k (x) in eq. (|5.21|) , we have 
X(y) = x, and using eq. (|5.8|) : 

E{x',Y k {x)) 



P(x,Y k (x);x') = d x 



i.e. 



P(x,Y k (x);x') = d x 



_(x - x')(Y k (x) -Y(x'))_ 
E(x',Y k (x)) 



+ 



-n(Y k (x);x')E x (x,Y k (x)) (5.22) 
E r (x, Y k (x)) 1 

(5.23) 



[(x - x')(Y k (x) - Y(x')) y{a k {x)){a k {x) - a{x'))_ 

We know that Pisa polynomial in y of degree d 2 — 1 , and we know its value in d 2 + 1 
points, therefore we can determine it with the interpolation formula (see appendix A), 
which, in this case, reduces to: 

P(r ... rp'\ _ V p ( x ^);x') E(x,y) 

h (V-Mx)) E y (x,Y k (x)) 

The latter expression seems to be a polynomial of degree d 2 in y, but one can check 
that the leading term vanishes (see appendix B). Thus we may also write: 

p( a = P(x,Y k (x);x')(Y k (x) -Y Q (x)) E(x,y) 
[,V,) h (y-Y k (x))(y-Y (x)) E y (x,Y k (x)) ^ 

Both expressions give either: 

d 2 



P(x, y; x') = d x , ^ — 
4 



E x (x,Y k (x))E(x,y)(Y k (x) -Y (x)) 



k J x y'Mx))(a k (x) - a(x'))E y (x, Y k (x))(y - Y k (x))(y - Y (x)) 
E(x',Y k (x)) E{x,y) Y k (x)-Y (x) 



(5.26) 



x - x')(Y k (x) - Y(x')) E y (x, Y k (x)) (y - Y k (x))(y - Y (x)) 



Using E(X(s),y(s)) = we have y'{a k {x))E y (x,Y k {x)) = -X'(a k {x))E x (x,Y k (x)), 
and we can write: 



P(x, y; x') = d x , -y, 



E(x,y) 



Yn(x) —Yu(x) 



< X'{a k {x)){a k {x) - a{x>)) (y - Y k (x))(y - Y (x)) 
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E(x',Y k (x)) E(x,y) Y k {x)-Y {x) 

(x - x')(Y k (x) - Y(x')) E y {x, Y k (x)) (y - Y k (x))(y - Y (x)) 

(5.27) 



or, using eq. Q5.24Q : 



P(x, y; x') = d x > J2 k=0 



d 2 E(x,y) 1 



X'(a k (x))(a k (x) - a(x')) (y - Y k (x)) 

E(x',Y k (x)) E(x,y) 1 

(x - x')(Y k (x) - Y(x')) E y (x, Y k (x)) (y - Y k (x)) 

(5.28) 

Eq.( |5.28| ) is clearly a polynomial in y, let us explain how one can check that it is 
also a polynomial in x: it is symmetric in the a k (x)'s, therefore it has no cut (it takes 
the same value in all sheets), it must be a rational function of x. The only possible 
poles could be at x = x', X\a k (x)) = or E y (x, Y k (x)) = 0. 

- Take x = x', more precisely, for some k e [0, d 2 ], cr k (x) = cr(x'), because of the 
symmetry one may assume that k = 0, and one sees that expression eq. ( |5.27| ) has no 
pole at ao(x) = a(x'). 

- Now choose some k and x such that X'(a k (x)) = e is very small. There must exist 
some (unique because our potentials are non critical) j ^ k such that (Tj(x) — a k (x) = 
0(e), and one has X'(aj(x)) = —X'(a k (x)) + 0(e 2 ). The poles for j and k thus cancel 
each other. 

- Now choose some k and x such that E y (x, Y k (x)) = e is very small. This means 
there must exist some (unique because our potentials are non critical) j ^ k such that 
Yj(x) — Y k (x) = 0(e), and one has E y (x,Yj(x)) = —E y (x,Y k (x)) + 0(e 2 ). The poles 
for j and k thus cancel each other in eq. (|5.28|). 



5.7 Determination of U (x, y; x') 



We can now compute U(x, y; x') using eq. ( |5.18|) : 
(y-Y (x))U(x,y;x>) = d x , £ — E{x ' v) 



- X'{a k {x)){a k {x) - a(x'))(y - Y k (x)) 
E(x,y)_^ E{x',Y k {x)) 



E V~* (5.29) 



(x — x')(y — Y(x')) 



which simplifies to: 



(y-Y (x))U(x,y;x>) = d x> ^— E{x ' v) 



- X'{a k {x)){a k {x) - a(x'))(y - Y k (x)) 
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E(x,y) _^ E(x',Y k (x)) 



E(x',y) - E(x,y) 



(5.30) 



(x — x')(y — F(x')) 
Note the following identity (see appendix B, Lemma |B.5| ): 

^ n? = o(n(x) - is-fro) _ ry*/-^)) r ,, n 

~ (Y k (x) - y) U^ k (Y k (x) - ^(y - Y}(x)) 1 ' j 

which implies: 

A £(x',n(x)) g(x^y) 

A.(y^)-y)^(x,y fe (x)) s(x,y) 1 ' j 

and thus: 



TT f a q ^(x,y) 
U{x,y;x) = d x i^ — 



- A"(<7 fc (x))Mx) - a(x'))(y - F fe (x))(y - F (x)) 



(5.33) 



It is easy to check that this expression is a polynomial in y and goes to when x —>■ oo 
in the physical sheet (sq — > oo) or when x' — > oo in the physical sheet. It has inverse 
square root singularities near the endpoints where X' (crfc(x)) = 0. 

Computation at y — Y (x) 

U( X , y oW ;*<) = 4,g x , {ak(x))(a ^f-J^ {Y<l{x) _ nW) < 5 - 34 > 

Computation at i' = x 

lX ' ° W ' XJ ~ ^ A"(a fc (x))^(a (x))(a fc (x) - a (x)) 2 (y (*) - F fc (x)) 
We are now equipped to compute the next to leading order functions... 

6 Next to leading order 

6.1 The one-loop function: computation of Y~W 

We have (see eq. (ggg)): 



(6.1) 
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where P^- 1 ' has degrees {d\ — 1, c?2 — 1), and the coefficient of x dl ~ 1 y d2 ~ 1 vanishes. 
We thus have d\d 2 — 1 unknown coefficients to determine. On the other hand, there 
are 2d\d 2 — 2 values of s for which E y (X(s),y(s)) = which are not endpoints (see 
appendix D, eq. ( |D.9| )). If we write that the poles at all these points cancel in equation 
|6.1| , we can determine P^(x,y). However, we don't need to determine P^(x,y), we 
will follow another method which allows to determine Y^'ix) directly. 

Eq.(|6.1|) is a rational function whose only poles are the endpoints, with degree up 
to 5. Since eq. ( |6.1| ) behaves as 0(s~ 2 ) when s — > oo, Y^ 1 ' has no pole at oo, and since 
eq. ( |6.1| ) behaves as 0(s l+d2 ) when s — > 0, X'(s)Y^ has no pole at 0. 

Therefore we may write: 

\r(i)(y( _\\ _ \ \^ Ak B k Ck Dk / fi 9 n 

[ * [S)) ~ X'(s) ^ (a - e k y + (s - e k f + (s - e k f + (s - e k ) ^ 

The coefficients A k , B k , C k , D k are determined by matching the poles in eq. ( |6.1| ) 
(P^ doesnot contribute to them). We will see below that D k = 0. This is not sur- 
prising since we expect 

yd) 

(x) to be a derivative (indeed Y^\x) = —d/dVi(x)F^ = 

-d/dx d/dV{(x)F<V). 

Let k G [0, d 2 ], and choose s close to e k : 

s = e k + e (6.3) 

there must exist s (unique because the potentials are non-critical) such that X(s) = 
X(s) and s is close to e k (s is the <J k (x) in eq. ( |6.1D ): 



s = e k -r] , 77 = 0(e) (6.4) 
By solving X(s) = X(s) order by order in e we get: 

r/ = Ae , \ = l+r k e + r 2 k e 2 + (2r 3 k + t k )e 3 + €>(€*) (6.5) 

where 

lX"'(e k ) 1 X IV (e k ) 1 X v (e k ) _ 

Tk ~Yx^ ' Sk ~ 6 x"(e k ) ' tk ~ eo;r(^y ~ TkSk (6 ' 6) 

From eq. ( |6.1|) we must have: 

Ak + Bk( + Ck<? + Dk , = _ _ _ _ _ _ _ _ _ _ + 0(e .) (6 . 7) 

We note the 3rd degree polynomial: 

P fc (e) = A fc + 5 fc e + C k t 2 + Ae 3 (6.8) 
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i.e. 



P fc (e) = (1 + A)- 2 

(-XX"(e k ) + Sx"\e k ) - e 2 ^X IV {e k ) + e^X v {e k )y l 
2 6 24 

((1 + A)y(e fc ) + 1(1 - \ 2 )y'\e k ) + f (1 + A 3 )y"(e fc ) 
2 6 

+^(l-A 4 )^(e fc ))- 1 + 0(e 4 ) 

= _(^( ejfc )y( Cjfc ))-i (l + A)-^- 1 

3 5 
(1 - -eXr k + e 2 X 2 s k - -e\t k + r k s k )y l 

d+ e (i \) r{ek) + e2 (i x + x 2 ) y "' {ek) )-i 

(1 + 2 (1 - A) 3^y + 6 (1 - A+A) 3^y } 

+0(e 4 ) (6.9) 

It is easy to see that t k as well as y IV (e k ) disappear, and we find that D k = 0. 
After a straightforward calculation, one finds: 

-8X"(e k )y'(e k )(A k + B k e + C k e 2 + D k e 3 ) 

{ _ 1 ,T"(e fc ) 1 2 ? A""(e fc ) 2 _ * JV (e fc ) *'"(e fc )y"(e fc ) _ y"(e fc ) \ 
3 e A"'(e fc ) 6 C \X"(e k f X"(e k ) X»(e k )y>(e k ) y(e k ) J 

(6.10) 

After substitution into eq. (S.2), we find the genus one correction to the resolvent: 



d2 -i 

1 K " f^8X'(s)X"(e k )y>(e k )(s-e k )± 

y> X"\e k ) 

f^ Q 24X>(s)X"2(e k )y>{e k )(s - e k f 

d2 X"'{e k f _ X IV (e k ) X"'{e k ) y"(e k ) _ y'"(e fc ) 
X"{e k y X"(e k ) f X"(e k ) y'(e k ) y'(e k ) 



E 



(6.11) 



, u A8X'(s)X"(e k )y(e k )(s-e k y 



That function represents the partition function of a statistical physics model on a genus 
one surface with one + boundary. 

6.2 The free energy 

We are now going to find the free energy F^ 1 ' such that: 

y<1) W - "av^) (6 ' 12) 
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For that purpose, we need to compute the derivatives of various quantities with respect 
to Vi(x). From our knowledge of the one matrix case |3|], we guess that F^ 1 ' will be 
related to the "moments" y'(e k ). Therefore, we shall compute the action of d/dV\(x) 
on such quantities. 



6.2.1 The derivatives of y(u) and X(u) 
In this section, we note x = X(s) and: 

X(u) = X'(s) dX{u) 



dVAx) 



y{u) = x (s) 

oVi[x) 



Notice that (using eq. O and = *™ = ™$Y' (*(«))+ §gg 



y\u)x(u) - x'{u)y{u) = —^—- 2 

In particular, if e^, k = 0, . . . , is an endpoint X' (e&) = we have: 

#(e fc ) ' 



y(e fc )(s-e fc ) 2 
Now, notice that A'(w) is a Laurent Polynomial^ in u: 



d 2 

X(u) = ju + a k u 

k=0 



where 



7 = X'{s) 



c?7 



a k = X'(s) 



da k 



k = 0, . . . , d 2 



dV^x) ' v 'dV^x 

We know its value at all e&'s, k = . . . d 2 , therefore (see appendix A): 

7 1 



X(u) = uX'(u) 



7 



fc=0 



e k (s - e k ) 2 (u - e k )X"(e k )y'(e k ) 



Since ad 2 = g* d2+ iy d2 , by comparing the coefficient of u 1 and it d '\ we find: 

d 2 



7 1 «rf 2 



1 



and thus: 



7 rf 2 « d2 2 ^ (s - e k ) 2 e 2 X"(e k )y(e k ) 



X'(s)^ = X(u) = ^Mf2 U + ek 



dV 1 (x) 



2 f u-e k (s- e k ) 2 e 2 X"(e k )y(e k ) 



(6.13) 

x'=X{u) 

(6.14) 
(6.15) 



(6.16) 



(6.17) 



(6.18) 



(6.19) 



(6.20) 



10 Indeed c? 2 is independent of V\(x). y{u) is not a Laurent polynomial, because the degree of y(u) 
depends on V\. 
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Notice that, as a function of s, this expression has double poles at the endpoints e k . 
For any i — 0, . . . , g?2, X'(ei) = implies that X'(s) d y^ = — *»(e\ anc ^ thus: 



X»( ei ) 2 f£e l -e k (s-e k yelX"(e k )y>(e k ) 



A"'(ei) 



2e i (s-e i ) 2 ^"(e i )y(e i 

From eq. ( |6.14j ) and eq. ( |6.20| ) we have: 



(6.21) 



which implies: 



aVi (x) 

X{e i )y'"{e i )+X"{e i )y\e: ° 



2AT"(ei)V V-V ■ " ( s _ e .)4 

+ ( S -e,) 3 A"'( ei ) J 

(6.23) 

This expression has poles of degree 2, 3 and 4, at the e^'s. We write: 



with: 

5k,i g _ &k,jX" (ej) 

'X"(ei) ' k,i ~ X"(e, 

and: 



Akl = -3-P^- , B ki = Z , V, (6.25) 



(^(e l )(y"(e l ) - y (e,)5g) - ^(ejyfo) jgffi + g'^y'je, y 
° Kl ~ 2Af"(e i )e^"(e fc )y(e fc ) (6 ' 26) 

where 

^'(m) m + e fc 
2 u — e k 



in particular: 



fcte) = e^"(e 4 ) , sK*) = h t X"( ei ) + \e\X"\e t ) (6.28) 



28 



and if k ^ i: 



g'He,) = X'\e t ) + ^'"(e.) + \e\X n '(e,) (6.29) 

9k{ei) = , 9' k (ei) = e i^ ± ^X"(e i ) (6.30) 

2 e 4 - e k 

d(ed = ^^X»(e t ) - 2 6 ^ k X"{ ei ) + e f-^X"\e % ) (6.31) 
We thus have: 

2*"(e i )y(e i ) 4A"'( ei ) 3 2e?*"( ei ) + e*"^) 2 1 ' ' 

and if k ^ i: 

(e 2 -e\- 2e iek )y> fe) 

2e|( ei - e k yX"(e k )y(e k ) ^ 66) 

Now, let us compute the following quantity (our guess is that it is related to Y^): 
y,, \ din yje-j) ^ A k , k ^ B k)k 



\ ^ ^k,k \ ^ \ ^ Ojy 

+ Z^ vrp,.v« - p,^2 + Z^ Z^ w 



(6.34) 

substituting eqs. (|6.25| , |6.32| , |6.33|) into eq. ( |6.34|) , we get: 
h dlny'( ei ) ^ 1 



X'(s) V ^ ^ - -3V 



,-_„ ^(*) f^X"(e k )y(e k )(s-e k y 

+ V_ *"'(e fc ) 



, X"{e k fy{e k ){s-e k f 

d 2 y'"(e k ) y"(e k )X'"(e k ) l X"'(e k f l X IV {e k ) 
, y(efc) y(e k ) *"(e fc ) 2 X"{e k Y J_ 3 X"(e k ) 

Z. 2^(e fc )y( eA; )( s - efe ) 2 



fc=0 

d 2 



2e 2 k X"(e k )y(e k )(s-e k y 

(ef - e 2 , - 2 ei e fc ) 



(6.35) 



t t i k 2el(e i -e k yX''(e k )y(e k )(s-e k y 



Notice the following identity^: 

5? < e -- e <=> -fe( e '-^ 2 



^ne proves it by taking the log derivative of u d2+1 X 1 (u) = Y\ k {u — e k ) 
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therefore: 

d 2 



dVy(x) 



We recognize: 



(6.36) 



2 ^(e fc ) e ^ w (e fc ) 2 



1 



" 3 £^)y(e fe )( S -e fe )4 

d 2 



*'"(e fc ) 



fc=0 



A"'(e fc ) 2 y(e fc )(s-e fc ) 3 



d 2 y'"(e k ) _ y"(e k ) X"'(e k ) _ X'"{e k f X IV (e k ) 

+ y(e k ) yin . ,v"(.» .\ •'(<,, i- ,\ " ( . - 



fc=0 



2X»(e k )y>(e k )(s-e k y 



(6.37) 



+ 1} S 2e^(e fc )y(e fc )( S -e fc ) 2 



(6.38) 



6.2.2 The free energy 

The free energy is such that: 

therefore to order 1/N 2 we have: 



<9F 



V{(x) - Y(x) 



av l{x) = - Y(1) M 

which implies (up to a constant K independent of Vi): 



where: 



/ —7 «! 




D 



1 



7 



dr+d 2 +2 



det 





diPdi 


V 







-7 



-7 
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\ 








A 









-7 / 



(6.39) 



(6.40) 



(6.41) 



(6.42) 



D is such that: 



fl=<-u* +i (^Y n ^^(-D-v n (^) 

^ ' / X'(e)=0 ' V'f?1=n \ ' / 



(6.43) 



y(g)=o 



The Constant K is independent of Vi, and if we repeat the calculation by exchanging 
the roles of M\ and M2, we find the same expression for F^- 1 ' with K independent of 
V 2 as well. By comparing with the gaussian case we get K = 0. 
We have found the genus one free energy: 



24 w y 



(6.44) 



Remark: D is a determinant which vanishes everytime the surface is singular, i.e. when 
the potentials are critical. We see that near a critical point, the genus one free energy 
has a logarithmic divergence. This is in agreement with known double-scaling limit 
exponents. 

6.3 Higher genus expansion 

The procedure can be repeated to higher orders. We have the three equations: 

E(x,Y(x)) = (V((x)-Y(x))(V 2 (Y(x))-x)-P(x,Y(x)) + l = -^U(x,Y(x);x) (6.45) 

N z 

(y - Y(x))U(x, y) = (V((x) - Y(x))(V 2 (y) - x) - P(x, y) + 1 - ±U(x, y; x) (6.46) 

Imagine that we already know the functions Y(x), E(x,y) (i.e. P(x,y)), U(x,y) and 
U(x, y; x') up to order h — 1. Let us write eq. (|6.45|) to order h. we know all terms but 
Y^(x) and P( h \x,y). The factor of the Y^ h \x) term is Ey°\x,Y^(x)) 1 it vanishes 
for all (x,y) such that y = Y^°\x) and Ey°\x,Y^°\x)), i.e. at 2d\d2 — 2 points, 
therefore we know P^ h \x,y) in 2did 2 — 2 points. Since P^ h \x,y) is a polynomial of 
degree (d\ — 1, ofo — 1), with diG? 2 — 1 unknown coefficients, P^(x, y) can be determined 
completely. Then, from eq. ( 6.45Q , Y^ h \x) can be determined. Using eq. ( |6.46| ), we 



can determine U^(x,y), and using eq. ( |6.47| ) we can determine U^(x,y;x') to order 
h. 

We clearly have an algorithmic recursive procedure to find Y^ h '(x) to any order h. 
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Finding the free energy is more difficult, as one needs to integrate with respect 
to the potentials. The easiest is to have an ansatz for and check that its derivative 
with respect to V\(x) is indeed -7^(4 as in ||. 

Define: 

M 0>i = y'fa) , M,.; = \ tJ k > 1 (6.48) 

Jo,i = X'\ei) , J k ,i = —, \ t} k > 1 (6.49) 

A. (e,J 

We expect the genus h > 1 free energy to be a rational function of these quantities, of 
the form: 

= E EE* m. •.^h nun ^ 1 j1 k '- iJ, '- < (^o) 

r,s,p,g,m>Oifci,...,fc 3 /i,...,/ i , ' i=0 °> l O' 4 

with the restrictions that: 

p + g = ^% + ^/ u (6.51) 

k u 

The selection rules (of the form m + p + q<4h — 4 as in ||) are still to be understood. 

7 The one matrix case: efe = 1 

In this section we check that for g?2 = 1, our results reduce to those of ||. 

Consider = 1 and V^y) = ^g^y 2 + S 1 *?/- ^2 is thus quadratic and M 2 can easily 
be integrated out. The 2- matrix model in this case reduces to a one matrix model with 
the potential: 

V(x) = VAx) - —x 2 + ^-x (7.1) 
The function X(s) and y(s) are: 

X(s)= 1S + a + ^ , y( s ) = ^ + Y j P k s k (7.2) 

s s 

k=0 

From eq. (|4.4|) we find that: 

«i = 92*1 , «o = fl£A> + S-i (7.3) 



V?(*(a)) = + - fa - 1 - (7.4) 

s 92 



and: 



indeed this expression is a polynomial in x and it satisfies eq. ( f4.4|) . It follows from 
eq. (O): 

V'(X(s)) = y(s) + y(g* 2 /s) - 2(/3 + I + -la) (7.5) 

s 9 2 
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The endpoints a and b are the roots of X'(s) i.e. they correspond to: 

e = ^/g* , a = X(-e) , b = X(+e) (7.6) 
The resolvent (defined in eq. ( |2.6| )) can be written [[|: 

W(x) = V((x) - Y(x) = - (v'(x) - M(x)y/(x-a)(x-bj S ) (7.7) 



where M(x) is a polynomial in x (the same as in 0). Using eq. (7J5) and eq. (|7.2|) we 
find that: 

M(X(s)) = y(S) ~ x ^ /8) , y/(X(s)-a)(X(s)-b) = sX'(s) (7.8) 
In particular at the endpoints we have: 

M(a) = ^ , M(i)=?M (7.9) 

7 7 



It is clear now, that our result eq. (|6.44 ) coincides with the genus one free energy 
found by the authors of 0: 

FW = -lln (^^y M(a)M(b)j (7.10) 

8 Conclusion and prospects 

The recursive method presented in this article, allows one to compute the resolvents 
as well as other loop-functions, to any order in the topological expansion of the 2- 
matrix model. We have found the expression of the genus one free energy, but an 
algorithmic method to compute the free energy to higher orders is still incomplete. For 
that purpose, it would be usefull to define some "moments" and basis functions, as 
in || , onto which one could decompose the loop-functions, and selection rules to know 
what set of moments should appear in the free energy at a given order. 

It would be interesting also to understand what terms survive in the double scaling 
limit when the potentials are near a critical point, and to develop a recursive method 
to compute directly the free energy and resolvents in the double scaling limit. 

As we said in the introduction, we have only developped a method to compute the 
terms in the 1 /N 2 expansion, but we have not proven the existence of that expansion. A 
rigorous proof of that existence can probably be obtained by Riemann-Hilbert technics 



It would also be usefull to understand what the 1/N 2 expansion becomes in the 
non genus-zero case. It is very likely that it would involve the ^-functions associated 



to the algebraic curve, as in [12 . 
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Appendix A The interpolation formula 

Let d be an integer, and Xq, . . . , be d + 1 distinct complex numbers. We note: 

d 

ir(x) = ]J(x ~x k ) (A.l) 

Let Po, . . . , Pd be d + 1 complex numbers. 

There is a unique polynomial P of degree d such that: 

P{x k )=P k k = 0,...,d (A.2) 

That polynomial is: 

d 

, (l {x - x k )7r'(x k ) 
Any polynomial P(x) of degree > d such that: 

P{x k )=P k k = 0,...,d (A.4) 

can be written: 



Tc(x) 

where Q(x) is a polynomial of degree degP — d — 1. 



^ {x - x k )7r'(x k ) 



Appendix B Proof that eq. ( |5.24j ) is a polynomial 

of degree di — 1 



The leading term of eq. (]5.24|) is (we use eq. (|5.28|) ): 



dx'^—Y 

nn /f> I ^^^^^ 



dl r /•;(./•'. V/.. (./•)) 1 x-x' 



X — X' 

k=0 



(Y k (x)-Y(x'))E y (x,Y k (x)) X'{a k {x)){a k {x)-a{x')) 
We have (see eq. ( p.5|) ): 



(B.l) 



d2 y 



iQ X'{a k {x)){a k {x)-a{x')) 



34 



Therefore we have to prove that: 



A E(x',Y k (x)) 1 
^ Q (Y k (x)-Y(x'))E y (x,Y k (x)) 



= 1 



(B.3) 



Note that: 



A E(x',Y k (x)) i _A n J - #0 (n(x)-is-« 
^ (n(x) - Y k ( X )) ^ rw 1 ^) - *;■(*)) 



(B.4) 



Lemma: 

For any Y , . . . , Y d2 and for any Y{, . . . , Y' d , 



r ' d we have: 



d 2 

E 



= 1 



(B.5) 



Proof: 



n^ol^fc _ i\k f(y s ^djyo, ■ ■ ■ , Yk-i, Yk+i, ■ ■ ■ , Yd 2 ) 

ttn^-Y,) ^ )T{k) " A d2+1 (Y , . . . ,Y d2 ) 



where f(Y) = Ylj^o(Y — Y-) is a monic polynomial of degree d- 
Vandermonde determinant of n variables. 



(B.6) 

! 2 in Y, and A n is the 



idermonde determinant ot n variables. 

We recognize a determinant (which is 1 by linearly combining columns): 

(\ Y Y 2 ... Y^- 1 

det 



d 2 

E 



1 Fo V 2 ••• Y d ^ f(Y )\ 
1 Y l Y, 2 ... Y^- 1 /(Yl) 

3 r d 2 • • • Y d 9 , 
Ad 2 +i(Y , • • • , Vdi) 



Vi ... iff" 1 /(Vcb)/ 



= i 



(B.7) 



Appendix C Other loop-functions 



Define: 



1 



H(x,y,x',y') = — ( tr 



1 



1 



2iV \ x - Mi y - M 2 x' - M x y' - M 2 
1 / 1 1 11 



+ 



tr 



2N \ y' - M 2 x' - M x y - M 2 x - M 1 




(C.l) 



35 



Notice that we have: 



H(x, y, x', y') = H(x', y', x, y) = H(x, y', x', y) 



(C.2) 



We also define: 

F(x,y,x',y') 

we have: 



2N 



x - Mi y - M 2 x' - M x y' - M 2 

l_/tr— 1 1 VW-VKMJ X 

2N\ y> - M 2 x' - M x y - M 2 x - M x / 



F(x,y,x',y') = F(x,y',x',y) 

Loop equations: 

The invariance under the change of variable 
8M l = 



(C.4) 



i i i i 

x-Mi y-M 2 x'-Mi y'-M 2 



+ h.c. implies: 



(y' - Y(x))H(x, y, x', y') = - W ^ y) W } X '> y \ l - W (x', y')) - F(x, y, x', y') (C.5) 

x — X 

The same equation with y <-> y' gives: 

1 



H(x,y,x',y 



(x - x') (y - y') 
E(x,y)E(x',y') 



E(x>,y)E(x,y>) 



(x'-X)(y-Y')(x-X')(y'-Y) 



(x-X)(y-Y)(x'-X')(y'-Y')_ 

E(x,y)E(x',y') 



(C.6) 



it follows: 

F(x,y,x',y') = 



(x - x')(y - y')(x - X)(x' - X')(y' - Y') 
E(x',y)E(x,y') 

(x - x'){y - y'){x - X'){x' - X)(y - Y') 
E(x',y)E(x',y') 

(x - x')(x' - X)(x> - X')(y - Y')(y' - Y') 



(C.7) 



Appendix D Some relations 



The following formula are usefull for many calculations in this article. They are easy 
to prove. 



E(x,y) = AY[{<Ti{x) - aj(y)) 



1,3 



with 



A = (-l) d2+ V+"/^ 1+ i = (-1)" 2+1 ^T 



(D.l) 



(D.2) 
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d.2 dl 

7 TT/ / \ \ ,,, \ Pdi 



X(s)-x=-^l[(s-a t (x)) , y( a )- y = ^-JJ( a -5- J -(y)) (D.3) 



s s 
i=o i=o 



(D.4) 



d,2 -, -. di . 

v - = V - 



i=0 



(D.5) 



£ w (s,y (*)) = -^^T73T II Wx) - &s(Yo(x))) 



A d2 



y^\x)X\a Q (x)) nil^W - (D.6) 

i^O j^O 



n(x) - Yotx) = " ^( y oW)) (D-7) 



Limits: 



ct 3 d2 3 a dl+1 

Ey {x{si y(s)) ^ -^s d ^ , e v {x{s), y{ 8 )) sZo -d2^^s^-^ 

(D.8) 

E y (X(s),y(s)) is a Laurent polynomial: 

E y (X(s), y(s)) = -^-P 2dld2 _ 2 ( s ) (D.9) 

where P2d 1 d 2 -2(s) is a polynomial of degree 2did 2 — 2. Therefore, E y (X{s), y{s)) = 
has 2did 2 — 2 zeroes. 
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